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The strong coupling constant a s {fi) is extracted from the fits to charmonium spectrum and fine 
structure splittings. The relativistic kinematics is taken into account and relativistic corrections 
are shown to increase the matrix elements defining spin effects up to 40%. The value of ct s ((J>) at 
low-energy scale /j, = f .0 ± 0.2 GeV was found to be a s (fJ,) = 0.38 ± 0.03(e:rp.) + 0.04(theor.) which 
is about 50% lower than standard perturbative two-loop approximation and is in good agreement 
with the freezing a s behaviour. 
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I. INTRODUCTION 

The strong coupling constant a s (Q 2 ) was measured in 
many experiments at large momenta where perturbative 
QCD can be successfully applied [Q. Due to these inves- 
tigations our knowledge of the fundamental constant A 
in QCD at high energy scale is now much more precise 
than before. At the low energy scale any extraction of 
a s from experiment is complicated by essential nonper- 
turbative (NP) effects and the r decay is assumed to be 
the only lowest-energy process from which the coupling 
constant a s can be cleanly extracted: for r decay NP are 
argued to effects give a small contribution to the hadronic 
ratio R T The averaged value of a s at the scale of r 
lepton mass is now a s (M T ) = 0.35 ± 0.03 [0. 

In the present paper we suggest another piece of low- 
energy data - very precise experimental measurements of 
spin splittings of Xc mesons. For reasons discussed below 
these data can give a unique information on the strong 
coupling constant a s {fi) at low energy scale, /i < 1 GeV. 
Even though NP effects are essential for IP charmonium 
state with the size as large as R « 0.65 fm the NP 
contributions will be shown to affect a s (p) in a simple 
and controlled way and the choice of NP parameters is 
strongly restricted by a fit to charmonium spectrum. 

We shall not consider here the hyperfine shift of Xc 
meson with regard to the center of gravity of 3 Pj multi- 
plet. This shift is small and negative, |A#j?| < 1 MeV, 
and small perturbative (P) and NP contributions cancel 
in A# p ||,Q . To explain this shift a very precise theory 
of NP interaction at small distances is needed which is 
absent by now. 

The experimental magnitude of tensor and spin-orbit 
splittings on the contrary are large enough (> 35 MeV) 
and known with the accuracy better than 2%. 

The spin structure of P states in heavy quarkonia was 
already investigated several times in the framework of 
potential approach 
a s (/j,) obtained in 



We quote here the values of 
] from the fit to charmonium fine 
structure: a s (1.22 GeV) = 0.386(> ^ m,m = 1.30 GeV) 
@ and a s \fi = m = 1.20 GeV) = 0.35 in §. Note 



that however different static interactions were used in 
HU , nevertheless the resulting a s (/i) values at the scale 
fj, « 1.2 GeV are rather close to each other and small in 
both cases. The calculated value of a s (1.2 GeV) is close 
to that of a s (M T ) at the energy scale that is significantly 
smaller than M T . 

In this paper an improved analysis of the fine struc- 
ture data will be presented. First, we take into account 
relativistic kinematics. The wave functions and all ma- 
trix elements will be calculated with the help of spinless 
Salpeter equation as compared to the nonrelativistic un- 
perturbed Hamiltonian used in |^||. For the charmo- 
nium IP state v 2 /c 2 ps 0.4 is not small and relativistic 
corrections are important. For example, whereas the size 
of IP state decreases only by ~ 8% in relativistic case, 
the matrix elements like (r~ 3 ), (r~ 3 In mr), denning spin 
splittings, increase as much as 40% for any set of param- 
eters. As a result relativistic corrections affect the value 
of a s (/z) extracted from experimental data. 

The second improvement is based on a simple obser- 
vation. We suggest instead of spin-orbit matrix element 
a = (Vso( r )) an d tensor matrix element c = (Vr(r)) to 
use their linear combination ry = 3/2c — a. This combi- 
nation has two remarkable properties: 

i) its perturbative part r\p contains only a 2 , radiative 
corrections since 0(a s ) terms cancel due to the relation 

4 1} =2/3aW; 

ii) the parameter r/p does not explicitly depend on the 
renormalization scale mass /i, whereas each of the second 
order terms, and ap', contains a term proportional 
to ln/i/m (see Sect. 2). 

Due to these features using rj instead of spin-orbit ma- 
trix element a is much more convenient while fitting to 
fine structure data. The experimental values of matrix 
elements a, c and rj can be easily calculated from Xc me- 
son masses |j| : 

a exv — 34.56 ± 0.19 MeV, 
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= 39.12 ±0.62 MeV, 



■q exp = 24.12 ± 1.12 MeV 



(1) 
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Note here that r] exp is relatively large (~ 24 MeV) with 
a small experimental error < 5%. The cited numbers, 
Eq. [l], slightly differ from the ones used in || because of 
recent change in Xco mass jjj . 

Our fitting procedure includes also a fit to charmonium 
spectrum. Here we prefer more refined fitting to mass 
level differences than to the absolute values M(nL) for 
the given state. The most important mass differences of 
levels lying below the open charm threshold are 



f = {m)- 1 |l.86723(r~ 3 ) - (r 3 lnmr)} , 



(8) 



M(2S) - M(IS) = 595.39 ± 1.91 MeV, 
M(IP) - M(IS) = 457.92 ± 1.0 MeV, 
Af (25) - M(1P) = 137.47 ± 1.77 MeV 

Here M(nL) is a spin-averaged mass. 



II. PERTURB ATIVE FINE STRUCTURE 
PARAMETERS 



(2) 



which has a weak dependence on parameters of static in- 
teraction and charm-quark mass. Practically in all cases 
f = 0.12 ± 0.01. 

The simple connection between a 2 (y) and rjp, Eq. (0), 
will be used later in our fit to the experimental fine struc- 
ture parameter rj exp , given by Eq. (fTl) . 



III. STATIC INTERACTION 

All spin effects in charmonium are very small as com- 
pared to the level masses or mass level differences in 
Eq. (^) and therefore can be considered as a perturbation. 
At this point the choice of unperturbed Hamiltonian is 
of importance. In j(| Schrodinger Eq. with the Cornell 
potential was used for this purpose. Here instead we take 
relativistic spinless Salpeter equation (SSE), 



Every matrix element a, c or r\ includes both P and NP 
contributions: c = cp + cnp, etc. First, P interaction will 
be discussed. The spin-dependent P interaction is now 
known in one-loop approximation only. Therefore our 
analysis can be only done with 0(a 2 s ) terms. In coordi- 
nate space the spin-dependent interaction including a 2 
corrections in the renormalization MS scheme was cal- 
culated in Q . Matrix elements of spin-orbit and tensor 
interactions for a number of flavors n/ = 3, can be easily 
found from the potentials 0] , 



ap 



,(2) (1) 



(2) 



2a 2 (n) 



j 1 
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- 3 Mn^- 



"p - 

1.58193(r- 3 ) + 2.5(r~ 3 ln(mr)), } 



c P 



Ji) _ 4 a s (n) 



. P +c p ,c P - 3 m2 {r 



4a 2 M 



3 7rm 2 

-3\ 



|4.5(r~ 3 ) 



In- 



up 

3.44916(r~ 3 ) + 1.5(r -3 ln(mr)) } 



(3) 



(4) 



(5) 



(6) 



In Eqs. a h matrix elements will be calculated for IP 
state. From Eq. (Q) and Eq. (Q) one can see that a,p and 

c P contain \nfi/m with the large coefficient whenever 
fx =/= m. However in the linear combination r/p — 3/2cp — 
ap these terms are cancelled and the following simple 
expression is obtained for rjp, 



_ 3 (2) (2) _ 2/q , 

VP - 7T c p - Op a sKN, 

2 nm 

where /o is the combination of matrix elements, 



(7) 



2^/p 2 +m 2 + V (r) ^ nl (r) = M nl ^ nl (r) 



(9) 



In the framework of potential model this equation with 
square root kinetic term was successfully used in many 
calculations of meson masses and properties during last 
twenty years ||[l(|. Recently it was deduced directly 
from QCD under a assumption of area law of Wilson loop 
in the framework of proper-time Feynman-Schwinger ap- 
proach |ll]. Therefore SSE, Eq. (||), cannot be consider 
as an ad hoc potential model, but rather on the same 
grounds as the QCD sum rules approach. 

In Eq. (0) the charm-quark mass m is a pole mass 
defined by the pole of quark propagator [7[|. T he pole 
mass m is related to the running mass in MS scheme, 
mini 2 ), as 



T'pole 



m(m 



4 a 5 (m 2 ) 

3 7T 



(K--) 



a s (m 2 ) 



(10) 



For c quark K w 14.0. From Eq. (|10|) one can estimate 
that for example for m = 1.4 GeV and a s (m) ~ 0.3 the 
relation m/m ~ 1.23, i.e. m is about 20 -f- 30 % larger 
than MS mass m. The static potential Vo(r) in Eq. (||) 
is taken here as Coulomb potential plus linear confining 
term, as it was done in 10], 



v o{r) = Var 

3 r 



(11) 



Here the vector coupling constant civ (a*) differs from the 
running constant a s (n) in MS scheme. The connection 
between them was found in |l2| and for nj = 3 it reads 



1.75 



.(/*)] 



(12) 
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In Eq. (|Tl]) cty(/i) will be taken at some fixed point /j, 
which in general can differ from the scale uq which de- 
fines spin splittings because Coulomb interaction behaves 
as whereas spin-dependent interaction behaves as 
r~ 3 and therefore is more sensitive to smaller distances. 
However, in our calculations it was found that the choice 
/iRi/ij together with the additional condition, Eq. (12), 



gives rise to a good description of both charmonium spec- 
trum and fine structure. 

In Eq. (|ll]) confining potential was chosen to be lin- 
ear at all distances, as it is done in most papers. This 
is in agreement with direct lattice calculations of static 
interaction |l4j], but is at variance with OPE and field 
correlator approach, which require that the NP interac- 

and tend to 



rfuonic 
I to be 



tion should start at small r as const ■ 
linear form or only at r > T g , where T g is the 
correlation length, measured on the lattice in |jl 
Tg^O.Sfm. 

Recently some arguments have been given in favour 
of additional linear potential at small r ||16fl and it was 
explicitly found from P-NP interference m |17|] to be of 
the magnitude which effectively resolves the discrepancy 
and confirm the linear NP potential at all r. 

Therefore we assume below that linear potential or is 
valid at all distances. To control the dependence of our 
results on the choice of string tension the parameter a 
will be varied in some interval compatible with a good 
description of mass level differences in charmonium. 

For pure linear potential NP contribution to spin-orbit 
potential is given by Thomas interaction for which 



o-np 



2m- 



■(r 1 ), a NP < 



(13) 



The value of cj^p is defined by T>\ -correlator [ [13|Jl8| 
which is small in lattice calculations |l5|], yielding the 
estimate < c NP < 0.3 MeV |l|]. This result is in 
a good agreement with direct lattice calculations of NP 
tensor potential (n]], where Vr(r) was found at the dis- 
tances 0.2 < r < 0.6 /to to be V T (r) < 1 MeV. Thus 
value of cmp is much less than \cinp\ ~ 15 4- 20 MeV 
and therefore can be neglected in JJnp- The theoretical 
error due to neglected term in Eq. [ll] is small and will be 
taken into account in our analysis. Then 



Vnp — 2° NP ~ aNP 



2m 2 



(14) 



To calculate wave functions of SSE the expansion of 
ij) n l in a series over Coulomb-type functions, suggested 
in ||, was used. The numerical calculations for IP 
state provided high accuracy (better than 10 -3 ) for ma- 



trix elements like 



}, etc.. Some of them 



are listed in Table | both for SSE and nonrelativistic 
Schrodinger equation for three sets of parameters. Set 
A with to = 1.2 GeV, a = 0.2 GeV 2 and a v = 0.35 
was taken from the paper @. From Table | one can see 
that the difference between relativistic (R) and nonrel- 
ativistic (NR) matrix elements is about 8% for square 



root radius, ~ 10% for matrix element (r 1 ), defining 
NP effects, Eq. (|l|), and very large, ~ 40%, for matrix 
element (r~ 3 ). This growth of (r~ 3 ) is due to the de- 
creasing IP-state size in relativistic case. 

One should also note here that mass difference 
M(2S) - M(1P) is large in NR case (~ 170 MeV), i.e. 
~ 20% larger than the experimental value 137 MeV , 
Eq. (||). When relativistic kinematic is taken into ac- 
count then excellent agreement with experiment can be 
easily reached for this mass difference, see Table ||. 



IV. FITTING CONDITIONS 

Our fitting procedure includes two conditions: (I) ?/ = 
r] eX p and (II) c = c exp . Using Eqs. (Q) and ( |l4| ) we write 
the first condition as 

"„(2) (2) , v _ 



V = 2 C p - " 



p ' 2m 2 v ' 
Vexp = (24.12 ± 1.12)10~ 3 GeV 



(15) 



With the help of Eqs. (||) and (||) this condition can be 
rewritten in another form, 



(u) ■ — = (24.12 ± 1.12)10~ 3 - 



2m 2 



A (16) 



From here a s (/tt) can be defined through the known num- 
bers, 



2/o 



(17) 



where m, a are fixed parameters and /o is given by the 
expression in Eq. (^]). Our calculations show that fo is 
almost constant, /o « 0.12 ± 0.01, independently on the 
choice of other parameters of static interaction, m, a and 
ay- 

As seen from Eq. (|l7|) the coupling constant is pro- 
portional to i/m and its value also strongly depends on 
\fK, where A in Eq. ([h]) is the difference between r\ exv 
and the absolute value of NP spin-orbit matrix element 
|ajvp|- As it follows from our calculations this differ- 
ence can become negative for some small quark masses 
(rn < 1.3 GeV) and large a w 0.2 GeV 2 whereas the 
l.h.s. of the Eq. (16) is always positive. Therefore the 
solutions with such small masses and large a must be 
excluded from our fit to the experimental data. For ex- 
ample, the Set A from [|| considered in the relativistic 
case yields the parameter \a^p\ = 26.53 MeV which is 
larger than n exp in Eq. ( TO) a nd hence inappropriate. 

With the use of Eq. ( |17| ) a s (fi) can be precisely de- 
termined for the given set of parameters to and a. At 
this stage the scale a still remains undefined but 0(a s ) 



terms, c'p and Op 1 ', which do not explicitly depend on 



can be calculated. 
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At the second step one can fit the second condition, since NP spin-orbit matrix element |ajvp|, proportional 



39.12±0.62 MeV, where c 



„W , J 2 ) 



'CNP- 



As it was discussed above NP term cnp is small, cjvp < 
0.3 MeV, and as compared to c can be neglected. With 



already calculated above, Eq. (0), the second condi- 



tion can be rewritten as 



„(2) 



4a s (/i) 



(r- 



(18) 



Using Eq. (o) it can be represented as 



^{( r -3) [3 .44«)U; 
6 irm z l 



4.5 In - 



1.5(r 3 ln(mr))} = 
(39.12 ±0.62)10~ 3 



4 a »M / -3v 



(19) 



All matrix elements and a s (/i) in Eq. (|19|) are already 



known for given m, a and ay , therefore Eq. ( 19 ) fixes /i at 



some value /io- It was found in our numerous calculations 
that the value /i does not coincide with the c quark 
mass m in general case, and only for very special choice 
of charm quark mass and static potential parameters ay 
and a the condition /i = m can be satisfied. That might 
be the reason why the choice /i = m apriori taken in the 
paper || has difficulties with simultaneous fitting of spin- 
orbit and tensor splittings (for IP state) with the same 
a s (/it) if c quark mass m = 1.8 GeV. In our fit, without 
the choice fi = m, when the Eq. (Jl|) and Eq. © are 
satisfied, we "automatically" get a = a exp . 

Hence for fixed m and a one finds a 8 (fi) and /i = /io 
satisfying experimental data. In the next Section we shall 
check whether this choice satisfies also data on the whole 
spectroscopy of charmonium. 



V. THE CHOICE OF PARAMETERS 

We shall discuss here only those fits with given m, a 
and ay which yield good description of spin-averaged 
spectrum in charmonium. As a result there appear some 
restrictions on the magnitude of a and ay, but the choice 
of pole mass m remains relatively arbitrary. Even if one 
takes MS mass, m = 1.30 ± 0.20 GeV , as it is commonly 
accepted |ij , and makes use of the Eq. ( JlOj ) , then the val- 
ues of the pole mass m can vary from m min w 1.2 GeV 
to m max w 1.8 GeV . Therefore it is important to impose 
additional physical restrictions on m. From our fitting 
procedure it is clear that for any mass m the correspond- 
ing a s (^o) and /io can be formally found. However the 
dependence of a s (p) on /i is different for different m. 
Some restrictions on the value of pole quark mass come 
from our fit to the fine structure data. 

All considered solutions with different m can be sepa- 
rated in three groups: 

i) Quark pole mass is small, m < 1.3 GeV. Then in 
Eq. (Ufl) the difference A 2 4 MeV is also small 



torn , is large. The value of A remains small even for 
the small a « 0.17 GeV. Then according to Eq. (|l|) 
a s (n) ~ VmA is small, a s (/i) < 0.20. As a consequence, 
0(a s ) terms, Cp and dp , are not large compared to 



the second order terms c 

„C0 > (1) 



p ' 

(2) 
P 



and 



(2) 



and in some cases 



c P ' > Cp' were obtained. For small m the calculations 
give large value of scale, /i > 3 GeV , and in many cases 
fi > 10 GeV, so that in all cases ln/x/m > 1. Note that 
the pole mass m < 1.3 GeV corresponds to the MS mass 
m< 1.15 GeV. 

Therefore for small m and any a the fine structure 
parameters a and c strongly depend on fi, and a\ cor- 
rections give large or even dominant contribution. Those 
solutions will be excluded in our analysis as unphysical. 

ii) The pole mass m is large, m > 1.6 GeV, which 
corresponds MS mass m > 1.3 GeV. Then in contrast 
to small m case the value Aw 10 -j- 12 MeV is large 
since the matrix element |ajvp| is smaller for large m. 
Therefore a s (fj,) ~ yfrnK grows large, ct s (fi) > 0.42 in 
our calculations. However, for large m the value of fi 
was found to be small, /j, < 0.7 GeV, so that | In fj,/m\ ~ 
1 is large again. As a consequence the negative term, 
proportional to ln/i/m in Cp\ Eq. (^), and a^p\ Eq. (Q), 
cancel positive contribution from other two terms, and 
as a result 0(a 2 s ) terms are numerically rather small. In 

some cases ay even becomes negative. 

Thus the second order terms have strong dependence 
on fi and one gets relatively large value of a s (fi) ~ 
0.4 -j- 0.5 at small mass scale /x < 0.7 GeV. Therefore 
in this case one can expect large contribution from a s 
corrections which are still unknown and hence any deci- 
sive conclusions about those solutions with large m are 
now impossible. 

Hi) The pole quark mass is in the range 1.4-5-1.56 GeV. 



For those masses /^-dependent term in c P and a}p' does 
not dominate and typically |ln/i/m| < 0.4 or 0.6 < 
/V TO ~ 1-0- J us t f° r such masses the best fit to char- 
monium spectrum was obtained. 

From two fits, to the spectrum and the fine structure 
splittings, we have found out the following restrictions on 
the choice of m and a, 

m = 1.48 ± 0.08 GeV, a = 0.178 ± 0.008 GeV 2 (20) 

For such m and a the Coulomb constant ay is "au- 
tomatically" fixed by the fit to charmonium spectrum: 
ay = 0.42 ± 0.04. Note also that the pole mass m in 
the range given by Eq. ( p0| ) corresponds to the MS mass 
m« 1.18 ±0.07 GeF. 

With the parameters from Eq. (^) the calculated val- 
ues of /io are in range, fiQ w 1.0 ± 0.2 GeV. At this 
scale, /io ~ 1 GeV, the extracted value of a s {no) (with 
experimental and theoretical errors) is 

a s (/io) = 0.38 ± 0.03{exp) ± 0M(theory) (21) 



(2) 
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0.75 < Mo < 1-2 GeV (22) 

Our calculations for two different sets of parameters 
are represented in Table ||. 

Set D with m = 1.4 GeV, a = 0.183 GeV 2 , and a v = 
0.39 was selected because for this variant mo as m was 
obtained, i.e. the term with ln/z/m does not contribute 

(2) (2) 

to Cp and a p . Also for Set D a good description of 
charmonium spectrum (see Table [n]) was obtained. From 
the fine structure fit the extracted a s {no) and /io values 
are 

a s (fi ) = 0.312, mo = 1-36 GeV as m (23) 

This a s (no) is small as compared with a s (M T ) extracted 
from t decay If we put the additional restriction: 

a s (m < M T ) > a s (M T ) = 0.35 ± 0.03 (24) 

then we have to conclude that for Set D the Eq. ( pi} ) is 
violated. 

To control the dependence of a s (fx) on NP parameter 
a we have calculated a s (/i) for smaller a, varying a till 
the description of charmonium spectrum was becoming 
poor. For example , for a = 0.17 GeV 2 M{2S) - M{IS) 
and M(1P) - M(IS) were already 30 MeV lower than 
its experimental values. For this a — 0.17 GeV 2 (m = 
1.40 GeV and dy = 0.39) the extracted a s (no) is in- 
creasing, 

a s (^ ) = 0.373, mo = 0.90 GeV, (25) 

these fitted values lie exactly in the range given by 
Eq. @ and Eq. ©. 

The best fit in our calculations was found for the set 
of parameters C, 

SetC: m= 1.48 GeV, 
a = 0.18 GeV 2 , a v = 0.42. (26) 

Note that the pole mass m = 1.48 GeV corresponds to 
m as 1.17 GeV. For Set C the excellent agreement for 
the mass level differences in charmonium (see Table ||) 
was obtained and 

a s ((i ) = 0.365 ± 0.027(ea;p), mo = 0.91 GeV (27) 

To estimate the theoretical error coming from NP effects 
we have varied a and analyzed the a s dependence on er. 
For the smaller a, a = 0.17 GeV 2 , a s ([i) is increasing, 

a 8 Oo) = 0.408 ± 0.027(exp), mo = 0.94 GeV (28) 

For this value of a the fit to spectrum is worse as com- 
pared to a — 0.18 GeV 2 , in particular mass level dif- 
ferences are about 5 -j- 10% less than the experimental 
ones, Eq. (Q). Therefore a s (fJ>) given by Eq. ( |28| ) can be 
considered as the upper limit for the coupling constant 
a s (fi) with m = 1.48 GeV. Averaging two numbers, 



Eq. (|27j) and Eq. ( |28| ) , yields the extracted value of a s (p) 
in Eq. ©. 

One should be reminded here that the nonrelativis- 
tic analysis in || discussed in Introduction has given 
a s (1.22) = 0.386 (and larger m = 1-22 GeV) close to our 
resulting value, Eq. J2l]), however, in relativistic case the 
fitted value of a s (/i) would be smaller for the same pa- 
rameters m, a and ay used in nonrelativistic approach 
m f§. Here we would like also to note a remarkable 
agreement between our result Eq. ( pl| ) and the value 
a s (Mo) = 0.38± 0.05(mo = 0.93 GeV) extracted from the 
best overall fit to 2P state in botomonium This 
coincidence is probably not occasional. Both systems, 
2P bb state and IP cc state, have exactly the same size: 

^{r 2 (bb)) 2P = v/(r 2 ( C c)) 1P = 0.63 ± 0.03 /to. 

The extracted value a s (0.9) ~ 0.38 turns out to be 
smaller then the corresponding value in PQCD. If we 
take for n/ = 5, A (5 ) = 237^4 MeV in two-loop 
approximation from recent compilation |l[ and calcu- 
late A^ 4 ) and A^ 3 ) from the matching conditions at MS 
quark masses, = 4.3 GeV and m c = 1.3 GeV, then 
A( 4 ) = 338±|? MeV and A< 3 ) = 384.4t 3 ^ MeV. Then in 
two-loop approximation a s (l GeV) — 0.54^^05 i s about 
a factor 1.5 larger than a s (l GeV) as 0.36±0.03±0.03 in 
our fit, the behaviour of perturbative a s (fi) is shown as a 
band in Fig. 1. Seven points from our fits with different c 
quark masses, 1.40 < m < 1.48 are also shown in Fig. 1. 

The variation of quark mass gives rise to changes in the 
values of \i and as a result we get the dependence of a s (m) 
on the scale m in low-energy region. It is important to 
underline that for IP state in charmonium the scale pa- 
rameter m lies in low-energy region, fi < 1 GeV , in all our 
fits (our seven points are shown in Fig. 1) a s (/j,) < 0.42 
is in agreement with Eq. (|T]). In all cases the extracted 
values of a s (ii) are below those in perturbation theory. 

How to interpret this result? We cannot accept the 
point of view that the relatively small value of a s (l GeV) 
found here points out that A^ 3 ) should be smaller than 
the cited above perturbative value of A^ 3 ) « 384MeV 
which corresponds to a s (M z ) = 0.119 ± 0.002. In our 
fit to the charmonium fine structure we have met a very 
specific case when a s (fi) is extracted at low-mass scale, 
M < 1 GeV. In this region the a s (/i) freezing phe- 
nomenon is expected |20[| which can affect the resulting 
values of the strong coupling constant, as is discussed 
below. 



VI. a s FREEZING 

Perturbative evolution of a s at large distances is phys- 
ically modified due to NP fields, which create confining 
strings (and hence NP mass parameter) and make an ef- 
fective large-distance cut-off in all loop integrals. This 
mechanism leads to a specific regime of ct s evolution at 
small Q, which is called freezing, and was studied both 
theoretically J2(J and in experimental analysis EJ. 
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Theoretically the freezing of the strong coupling con- 
stant can be deduced if one takes into account the be- 
haviour of running coupling constant in background vac- 
uum fields |2(it| . Experimentally phenomenon of freezing 
was anticipated already for some time (see review in [ pl| ) . 

The coupling constant, denoted as as{q 2 ), appears to 
be an analytical function at all Euclidean momenta, be- 



cause it depends on q 



'B ' 



where m 2 B is a background 



mass. At large q 2 3> as(q 2 ) coincides with perturba- 
tive a s (q 2 ). It means that A^ = A^ 5 ) and as{q 2 ) keeps 
the property of asymptotic freedom. At small q 2 , as{q 2 ) 
is freezing at some constant value, which depends on the 
background mass tub [p2[ . The behaviour of as(q 2 ) in 
two-loop approximation is given by the following approx- 
imate expression |20 22]j2^ |: 



UB{q 2 



47T 



fh In t E 
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(29) 



with tB = \u^^ 

The background mass raj in Eq. (E9f) can depend 
on the process considered and for static interquark in- 
teraction tub coincides with the lowest hybrid mass 
|^0|,^2| . Analytic and lattice calculations jMj predict 
the energy of the lowest hybrid excitation in the inter- 
val 1 -T- 1.5 GeV . In that follows we find the best fitting 
value tub = 1.1 GeV, which is well inside the predicted 
interval. 



Since Av = A (5) and m B is fixed, A^ and A 



(3) 



can 



be calculated from the matching conditions at the fla- 
vor thresholds (nTf, = 4.3 GeV and m c = 1.3 GeV). 



The obtained values of A^ and A^ ; turn out to be very 
close to perturbative values. In two-loop approximation 
A^ } = 339±lf MeV is only by 1 MeV larger than A' 



(3) 



(4) 
B 

and A^ 3) = 400±|f MeV is by 16 MeV larger than A^. 

The freezing value of ctn(O) = O.S^q'qs can be found 
from the definition Eq. ( p9|) with tub = 1.1 GeV. The 
behaviour of as(p) as a function of \i is shown in Fig- 
ure 1. by the dashed curve. From Figure 1 one can see 
that all seven values of a s (/io) found in our fit with differ- 
ent m lie on this curve ag(/i) in the region /i < 1.2 GeV. 

At the scale of the r lepton mass the value of 
a B (M T ) = 0.306t°;°i4 is close to the lower limit of 
a s (M T ) = 0.35 ± 0.03 from the compilation Q. How- 
ever, in some theoretical approaches, e.g. in renormalon 
chain model, the value of a s (M T ) — 0.305 is preferred 
||], which coincides with our prediction. 

Note also that a s freezing phenomena is actually ob- 
served in recent lattice calculations of three gluon vertex 
function in where a s (fi) points at /i < 1.8 GeV lie 
far below than perturbative ct s (fi) values and all points 
look like going to some constant value. 



VII. CONCLUSION 

Two fits - to mass level differences and the fine struc- 
ture splittings in charmonium were done here. In all 
calculations the relativistic kinematics was taken into ac- 
count which is more important for spin-dependent effects 
than for the spin-averaged cc spectrum. 

The fine structure parameters experimentally known 
from precise measurements of \c masses are not small (~ 
35 MeV) and measured with accuracy better 2%. The 
value of linear combination rj = 3 /2c — a is also not small 
and known with 5% accuracy, the use of this combination 
allows to simplify an extraction of a s (/x) from the fine 
structure data. 

For the fixed charm-quark mass m and NP static 
interaction the value of a s (^o) an d the scale mass hq 
are unambiguously calculated from our fits. In practice 
the variation of string tension is admitted in the range 
0.17 < cr < 0.185 GeV 2 till a good fit to mass level differ- 
ences is obtained. The uncertainty in a value gives rise 
to the theoretical error ~ 5% in a s (fi). 

The dominant theoretical uncertainty comes from the 
choice of c quark pole mass which has broader range than 
for MS mass, m(m 2 ). Therefore we had to solve spin- 
less Salpeter equation for many sets of parameters with 
different m. 

The value of m = 1.48 ±0.08 GeV was obtained in our 
best fit which corresponds to the MS mass fn = 1.18 ± 
0.07 GeV and the extracted value of the strong coupling 
constant is a s (po) = 0.38 ± 0.03(exp) ± 0.04(theor) at 
scale /io = 1.0 ± 0.2 GeV The uncertainty connected 
with NP effects is included in theoretical error. This 
value of a s (l GeV) is rather close to that calculated in 
bottomonium for 2P state Q], probably, because 2P bb 
state has the same size as IP state in charmonium. 

We suggest to consider relatively small value of 
a s (l GeV) as a manifestation of a s freezing and com- 
pare our calculated values of a s (/i) at fj, ~ 0.7 1.2 GeV 
with theoretical formula of as (/x) in background pertur- 
bation theory. 

The behaviour of «s(g 2 ) at large energy scale coincides 

with perturbative predictions since A^ 5) = A( 5 ) that gives 
rise to a s (Mz) — 0.119 ± 0.002 in two-loop approxima- 
tion. For n f = 3 A^ = 400±^ MeV is slightly larger 
that perturbative A^ 3 ) = 384.4 MeV. With this 
the predicted value of a s (M T ) w 0.306^q;q24 at the scale 
of t lepton mass is smaller that the conventional value 
0.35 ± 0.03 but coincides with the prediction in renor- 
malon chain model. 
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FIG. 1. The running coupling constant as a func- 
tion of the scale parameter fi. The perturbative two-loop 
a a (ji) with A (5) = 237124 MeV , A<4) = 338±^ MeV and 
A (3) = 384.4+32 Me y ig ghown by a hatched band. The 
running two-loop coupling constant in background pertur- 
bation theory with the background mass raj = 1.1 GeV, 
A^ 5) = A (5) , A (4) = 339.2±|f MeV and A (3) = 399.61^ MeV 
is shown by a dashed lines. The values of a s {fio) extracted 
from the fit to charmonium fine structure data for different 
quark masses are depicted by open circles with error bars in- 
cluding both theoretical and experimental uncertanties. 
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TABLE I. The comparison of IP state matrix elements for spinless Salpeter and Schrodinger 
Equations. 



matrix 
element 


Set A 
m = 1.20 GeV 
a = 0.20 GeV 2 
a v = 0.350 


Set B 
m = 1.40 GeV 
a = 0.17 GeV 2 
& v = 0.373 


Set C 
m = 1.48 GeV 
a = 0.18 GeV 2 
& v = 0.420 


R a) 


NR a) 


R 


NR 


R 


NR 


{r- 3 ) 1P (GeV 3 ) 


0.117 


0.080 


0.118 


0.086 


0.1423 


0.1014 


(r- l ) 1P (GeV) 


0.382 


0.345 


0.382 


0.351 


0.4040 


0.3704 


aNp(MeV) 


-26.53 


-23.96 


-16.57 


-15.22 


-16.60 


-15.22 


y/^ipiGeV- 1 ) 


3.281 


3.565 


3.293 


3.514 


3.130 


3.340 



a) R (NR) refers correspondingly to relativistic (nonrerativistic) case. 



TABLE II. The spin-orbit and tensor IP state splittings and spin-averaged mass level differences 
in charmonium a ' (in MeV). 





Experimental 
values > 
(MeV) 


Set C 
m = 1.48 GeV, 
o = 0.18 GeV 2 , 
a s (p Q ) = 0.365, 
Ho = 0.909 GeV, 
a v = 0.42 


Set D 
m = 1.40 GeV, 
a = 0.183 GeV 2 , 
a s (no) =0.312, 
Ho = 1-357 GeV, 
av = 0.39 


ClNP 




-16.60 


-18.40 






47.56 


41.84 






3.61 


11.12 


O-total 


34.56 ± 0.19 


34.56 


34.56 


c (1) 




31.70 


27.89 


c (2) 




7.42 


11.23 


Ctotal 


39.12 ±0.62 


39.12 


39.12 


M(2S) - M(1S) 


595.39 ±1.91 


591.30 


586.09 


M(1P) - M(1S) 


457.92 ± 1.0 


460.68 


445.93 


M(2S) - M(1P) 


137.47 ± 1.77 


130.62 


140.15 



a) all matrix elements denning spin effects were calculated for spinless Salpeter equation. 

b) M(nL) means the spin-averaged mass. 
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